Anomalous Quantum Diffusion at the Superfluid-Insulator Transition 
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We consider the problem of the superconductor-insulator transition in the presence of disorder, 
assuming that the fermionic degrees of freedom can be ignored so that the problem reduces to one 
of Cooper pair localization. Weak disorder drives the critical behavior away from the pure critical 
point, initially towards a diffusive fixed point. We consider the effects of Coulomb interactions and 
quantum interference at this diffusive fixed point. Coulomb interactions enhance the conductivity, 
in contrast to the situation for fermions, essentially because the exchange interaction is opposite in 
sign. The interaction-driven enhancement of the conductivity is larger than the weak-localization 
suppression, so the system scales to a perfect conductor. Thus, it is a consistent possibility for 
the critical resistivity at the superconductor-insulator transition to be zero, but this value is only 
approached logarithmically. We determine the values of the critical exponents r],z,v and comment 
on possible implications for the interpretation of experiments. 

PACS numbers: 74.76.-w, 74.40.-|-k 



I. INTRODUCTION. 

In a perfectly clean system at T = 0, the free Fermi gas 
is perched precariously at a critical point. An arbitrarily 
weak interaction will drive the system superconducting 
(by the Kohn-Luttinger effect if the interaction is repul- 
sive). In the presence of disorder, however, the diffusive 
Fermi liquid is a stable phase for a finite range of in- 
teraction and disorder strengths in dimensions d > 2. 
In d = 2, it remains an open problem whether or not 
fermions have a stable diffusive metallic phase. Such a 
phase, if it exists, cpuld not be adiabatically connected 
to the Fermi liquida since the non-interacting Fermi gas 
is always insulating in the presence of disorder in d = 2. 
In the limit of weak disorder, this can be understood as 
a quantum interference effect which is singular as a re- 
sult of the diffusive nature of electron propagation in a 
disordered system: diffusion at intermediate length scales 
(longer than the elastic mean-free path) thwarts diffusion 
at long scales (longer than the localization length)a. The 
interacting-electron problem remains unresolved because 
interactions in the spin-triplet channel are also si: 
as a result of the languid pace of diffusive motio 
The upshot of the interplay between these different sin- 
gularities is unknown (see, however, rcf. |^,^. 

Consider the critical point separating the insulating 
and superfluid phases of a perfectly clean system of 
bosons at r = in 2D. We would like to draw an analogy 
between it and the free Fermi gas. In the bosonic case, 
there is a particular value of the chemical potential for 
which the system has gapless critical modes, loosely anal- 
ogous to the excitations of the free Fermi gas. For any 
other value of the chemical potential, the bosons are ei- 
ther in a superfluid state - a superconducting state, if we 
assume that the bosons are Cooper pairs - or in a gapped 
insulating state. Suppose we now add disorder to this 
system. What is the fate of this critical point? On gen- 
eral grounds, we believe that it is unlikely to broaden into 



a stable diffusive metallic phase, and that the only stable 
phases are insulating (Mott insulator or Bose glasaj) or 
superconducting. Instead, we expect a diffusive metal- 
lic critical point with a universal conductivity separating 
the insulating and superconducting phases. The analogy 
between Fermi and Bose systems is imprecise, but it em- 
phasizes the important point that in both cases there is a 
ballistic critical point in the clean system which must be 
usurped by a diffusive fixed point in the disordered one. 

Such a fixed point should be amenable to analysis by 
methods similar to those used for the diffusive Fermi liq- 
uid. Conversely, expansion about the pure critical point 
- which is ballistic, not diffusive - should fail. In consid- 
ering such a perspective, one is faced with the following 
question: why do quantum interference effects, which ap- 
pear to be such an inevitable consequence of diffusive mo- 
tion, not preclude a finite conductivity at the superfluid- 
insulator transition? The answer must lie in the effects 
of interactions, which one might hope to tame since spin- 
less bosons, such as singlet Cooper pairs, do not have a 
triplet channel - the troublesome, singular one - through 
which to interact. 

In this paper, we present the results of such an anal- 
ysis. We find that there are two competing effects at a 
putative 2D diffusive Bose liquid critical point: one re- 
sulting from interactions between the bosons; the other, 
from quantum interference, i.e. weak localization. In 
the fermionic case, it is advisable to consider quantum 
interference and interactions on the same footing since 
they lead to similar logarithmic corrections at the per- 
turbative level. In the bosonic case, one must perforce 
do so, since quantum interference leads to the existence 
of localized states even in the weak disorder limit, and 
bosons would congregate in the lowest energy localized 
state in the absence of interactions. We find that the ef- 
fect of interactions is stronger than quantum interference 
and drives the system to a perfect conductor, thereby 
explaining how diffusion can remain impervious to local- 
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ization. This result is congenial to one's intuition that 
repulsive interactions should disfavor localization. Po- 
tential wells due to impurities diminish in attractiveness 
when they are occupied and, as a result, the random 
potential is effectively screened. This effect is present 
for both short-ranged interactions as well as long-ranged 
Coulomb interactions, but is stronger in the latter case. 
The same phenomenon occurs in fermionic systems as 
well, but it competes with the exchange part of the inter- 
action, which is opposite in sign due to Fermi statistics. 
If the interaction is short-ranged, it is irrelevant for spin- 
less fermions, so it has no effect on the conductivity in 
the infrared limit. (This is clear in the (5-function limit, 
where the direct and exchange interactions cancel.) In 
the case of Coulomb interactions, the exchange interac- 
tion between spinless fermions dominates and suppresses 
the conductivity. In the case of spin-1/2 fermions, the 
runaway flow of the triplet interaction amplitude indi- 
cates that the Hartrce interaction begins to prevail over 
the exchange interaction at longer length scales, thereby 
leading to an enhanced conductivity. However, the in- 
teraction strength diverges before a metallic fixed point 
is reached, and no conclusion can be drawn about the 
existence of a metallic state at zero-temperature. These 
difficulties do not arise in the bosonic case. The exchange 
interaction has the same sign as the direct one, and both 
enhance the conductivity. 

Our result is valid for large conductivities in units of 
/h. Hence, if the bare conductivity is large - as it 
can be if the bosons have an anisotropic mass tensor - 
then the renormalized conductivity is infinite. If the bare 
conductivity is small, then there are two possibilities. If 
the conductivity initially flows to sufficiently large values 
that we can apply our calculation, then it will continue 
to flow to infinity. However, it is also possible that the 
system will flow in this case to a different fixed point 
at which the conductivity if finite. In such a scenario, 
there would be two different possible universality classes 
of superconductor-insulator transitions. In either case, 
we conclude that it is a consistent possibility for the crit- 
ical point between the superfluid and insulating states of 
a disordered Bose liquid to be a perfect conductor. 

We derive these results in a non-linear cr-model 
(NLaM) formulation of the problem of diffusing, inter- 
acting bosons. Our NLoM is very similar to Finkel- 
stein's model for fermionscl. However, the NLcrM plays a 
very different role in this problem than in the fermionic 
problem. There, the NLaM describes the entire metallic 
phase. In 2-|-e dimensions, the metal-insulator transition 
occurs near the metallic fixed point, so the NLctM exconi- 
passes it as well. In the bosonic problem which models 
the superconductor-insulator transition, our NLcrM de- 
scribes the critical point. The antiferromagnetic Heisen- 
berg model in d > 2 provides an enlightening analogy. 
For isotropic exchange coupling = Jx,y, the model 
is ordered and is described by a NLcrM. In the ordered 
phase, continuous symmetries are broken so there are 
Goldstone modes; this is the analog of our critical point. 



For Jz > Jx,y, the model develops Ising order with a gap; 
this is analogous to our insulating phase. For Jz < Jx,y, 
the model develops XY order, which is analogous to our 
superconducting phase. 



II. DIRTY BOSONS 

Following ^, we will treat the Cooper pairs in a dirty 
superconductor as bosons moving in a random potential. 
We will assume that all fermionic degrees of freedom are 
gapped or localized and are therefore unimportan t .m T his 
assumption has been called into question recentlyll3llil. If 
fermionic degrees of freedom prove to play an important 
role at the superconductor-insulator transition, then our 
analysis will need to be modified to include them, but 
our description of dirty bosons will remain an important 
component of a richer description of the superconductor- 
insulator transition. r-. 

Note that we are studying here the generic transitions 
between the Bose Glass and superfluid phases which oc- 
curs at an incommmenusurate boson density. In the 
special case in which there are an integer number of 
bosons per lattice site, there may be a direct transition 
between Mott Insulating and superfluid phases which is 
tuned by varying the ratio of the hopping and interaction 
parametersE3. 

We begin with a system of interacting bosons moving 
in a random potential in two dimensions. The derivation 
which follows goes through in arbitrary dimension with 
minor changes, but d = 2 is the most interesting case. 
The imaginary-time action is: 

j (fx dT tl;* (^dr - ^ - ^ + V{x)^ V 
+ [ d^xd'^x' dTiP*{x)iP{x)u{x-x')iP*{x')ijj{x') (1) 



u{x — x') is the interaction between bosons; we will con- 
sider the cases of both short-ranged interactions and 
Coulomb interactions. V(x) is the random potential; we 
use the replica trick to average over it, thereby obtaining 
the action: 

Jd^x dr rai^, t) (^dr - - r) 

d^xdrdr' ^vq tp*{x,T)^pa{x,T)il;l{x,T)^pt{x,T') 
+ I d^xd^x'dT'iPl{x)^paix)u{x~x')il:*{x')ipa{x') (2) 



a = l,2,... , TV is a replica index. We have assumed that 
the potent ial has the Gaussian white-noise distribution 
V{x)V{x')^voS{x-x'). 

This action is problematic because it is not positive 
definite as a result of the second term. To cure this, we 
will rotate the integration contour in the functional in- 
tegral, as one does in the non-interacting case. This can 
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be done more conveniently if we work in the Matsubara tive definite by rotating the fields in the following way: 

frequency representation and separate the real and imag- ipnaA e~5**'sn(ri)0^^^^ A = 1,2. We rotate the n = 

inary parts of the Matsubara fields ipna = 4>nai + i4'na2, mode along with the n > modes, 
where e„ = 27rn//3. The action can be made posi- 

I 



The action now takes the form: 



+ ^ J -Vq (l)naAAnn'4>n'aA(l>mbBAmm'4>' 

n ti' ryi -m ' 



+ J2 d^xd?x 



^-iriY, sgn(mi)/4 



{x) u{x - x')(j),niaB{x')(j)miaB{x') (3) 



where 



sgn(m) Sr, 



r 



In the absence of disorder, repulsive interactions are 
marginally irrelevant, and the critical behavior of (||) is 
controlled by the Gaussian fixed pointla. Now consider 
a perturbative treatment of the disorder. In the self- 
consistent Born approximation, we find a self-energy due 
to disorder of the form: 



S(e„) 



mvo 



In 



AV2m 



iHn + /i + S (e„) 



i tan 



e„ + ImS (£„) 
/i + ReE (e„) 



(4) 



The random potential shifts the chemical potential and 
also gives the bosons a finite lifetime r. As a result of 
the lifetime r, single-boson excitations are no longer long- 
lived degrees of freedom. However, particle-hole pairs are 
long-lived, as may be seen from the conductivity which, 
at this level of approximation, is cr = . 

This does not preclude critical behavior in the single- 
particle properties, as has already been seen in the 
context of interacting fermionsEI and of quasiparticles 
in a disordered d-wave superconductor where there are 
density-of-states corrections and also in the context of 
non-interacting electrons with an extra sublattice sym- 
metry, pWjhere the single-particle Green function itself is 
criticallla. 

The conductivity is small because there are no particle- 
hole pairs for t = oo (since the transition occurs at 
the bottom of a quadratic band). A finite lifetime 
leads to a small density of states ~ 1/r for particle- 
hole pairs, which cancels the factor of the lifetime to 
which fj is customarily proportional, thereby leading to 
a conductivity which is 0(1). However, we note that a 
parametrically large conductivity can be obtained in a 
slight generalization to a model of two species of bosons 
with anisotropic masses and that mix upon scattering. 
Suppose that one of them has rux = mi, my = m2, 
while the other has masses reversed. Then we find that 



a — \Jm\lm2 4- \/m2lm\ /27r^. For sufficiently large 

or small ratio mi/rn2, the conductivity will be large. 
Such a situation could occur, for instance, in a two- 
band model in which the two bands of electrons have 
anisotropic masses, leading to anisotropic masses for the 
Cooper pairs. 

An RG analysis of the dirty boson problem yields the 
following RG equation in an e-expansion about d — 4lI: 



d£ 



—— = (e + er) + Bvfj 



(5) 



with 4 — e — Er spatial dimensions and Et time dimensions 
(the interesting case d = 2 occurs at e = e,- = 1)- B > 0, 
so there is no fixed point at weak coupling; instead, there 
is a runaway fiow to strong disorder. We interpret this as 
an instability of the pure critical point, at which the criti- 
cal modes are ballistic, to the diffusive fixed point. To ac- 
cess the latter fixed point, we will construct a non-linear 
a model which is appropriate for physics at length scales 
longer than the mean-free path. In this regime, transport 
is diffusive, and we may neglect degrees of freedom, such 
as the (j) fields, which are short-lived. 



III. SADDLE-POINTS FOR DIRTY BOSONS 

In the absence of the ie„ term, the non-interacting 
part of the action (||) has an O ((fc -I- l)iV, kN) symmetry, 
where fc is a cutoff on the Matsubara frequencies. The 
key assumption of Finkelstein's theoryEl for fermions is 
that the elevation of the energies of the diffusion modes 
by the ie„ term and the interactions can be neglected 
compared to the gaps associated with other degrees of 
freedom; when this condition is satisfied, it is valid to re- 
tain only interacting diffusion modes and ignore all other 
degrees of freedom. We make the same assumption here 
in our description of the critical point. In the superfluid 
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state, this is clearly not sufficient, and we will have to re- 
tain an extra degree of freedom. It may also be necessary 
to include extra degrees of freedom to properly describe 
the Bose glass insulating state. 

Our treatment of the critical saddle-point and non- 
linear CT-model (NLcrM) for interacting bosons follows 
that of Finkelstein for the fermionic case and also that of 
the bosonic representation of the non-interacting prob- 
lem. Hence, we will merely give an outline in this section 
and the next, emphasizing the important differences. De- 
tails are presented in appendix W. 



We begin by using the Hubbard- Stratonovich transfor- 
mation to decouple the Vq term with a matrix QJJ^"^^. 
We then decouple the interaction in two different ways 
with X, which decouples the direct and exchange chan- 
nels according to X ^ ^nd Xc, which decouples the 
Cooper channel according to Xc ~ V'V'- Finally, we de- 
couple the chemical potential term with $ ~ -0. In this 
way, we have a system of non-interacting bosons at zero 
chemical potential - their critical point - moving in the 
background fields X, X^ and Integrating out the (j) 
fields, we obtain the eS'ective action (see appendix ^): 



S,ff[Q,Y,Z, = 

E 



tr In ( ien + — 
2m 



2vo 



(6) 
(7) 

(8) 

(9) 



r 



The Green function G of the is written as a 2 x 2 
matrix in the final line to emphasize the particle-hole 
structure. It is the operator inverse of the expression 
inside the logarithm. For /i < 0, it is not even necessary 
to introduce $; we can simply drop the last line of (||) 
and insert fi inside the logarithm. 

Let us now consider the saddle-points of this effective 
action. For fi > 0, there is a saddle-point with ($) ^ 0. 
(When we include fluctuations, fj, will be renormalized, 
so the critical value will not be zero.) When $ develops 
an expectation value, Q, Xc, and X are forced to follow 
since they are coupled directly to bilinears in <&. This is 
the superfluid phase. 

For /i < 0, let us consider the non-interacting case 
r = Fc = 0. The saddle-point condition is 



Q 



1 



(10) 



Let us absorb the real part of the saddle-point value of Q 
into a renormalized /jr and focus on the imaginary part. 
For fiR — 0, the saddle-point solution of ( px| ) is 



b,AB 



« — sgn(e„) d„ 



Sab Sab 



(11) 



This is the diffusive saddle-point for self-consistent Born 
scattering of critical bosons by impurities. It corresponds 
to a finite density of states for the bosons at this level of 
approximation. Notice that this saddle-point solution is 
taken to be replica symmetric. 

Now, for < 0, there is another translationally- 
invariant saddle-point with Q = 0. For this solution, a 



non-zero density-of-states in not generated in the insulat- 
ing state at this level of approximation; it remains a Mott 
insulator. We would like to point out two possible mech- 
anisms to generate the finite density of states that occurs 
in the Bose-glass phase. One is that the correct saddle 
points are replica symmetry broken mixtures of the Q = 
and Eq. |ll| solutions. A possible self-consistent solution 
is one still diagonal in replicas but with zero matrix ele- 
ments for p replicas and unit matrix elements for n — p 
replicas. Another possibility is that there are non-trivial 
instantou-saddle-points which generate a finite density- 
of-statesEj. In the absence of interactions, the bosons will 
condense into these localized states, so we must consider 
the corresponding instantons with r,rc ^ 0. At present, 
we do not have a description of the Bose glass insulator, 
but this does not affect our ability to describe the critical 
point between it and a superconductor. 

It is useful, in thinking about this theory, to imagine 
lowering the temperature of a system of dirty bosons. 
At finite temperature, there will be a finite aredge in the 
phase diagram - the quantum critical regiontilEj - where 
the bosons will be effectively critical. In this regime, we 
may begin by considering non-interacting bosons which 
are semiclassically scattered by impurities. As we de- 
crease the temperature, we must begin to include the 
effects of interactions and of quantum interference pro- 
cesses. If we stray too far from the critical fi as we lower 
the temperature, thereby leaving the quantum critical re- 
gion, then we cannot include these effects perturbatively. 
It is clear that they completely destabilize the diffusive 
saddle point, so they must be included right from the 
start (e.g. by starting from new saddle-points, as we 
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have sketched above) in order to describe the superfluid 
or insulating phases correctly. However, so long as we 
remain at criticality, we can hope to account for these 
effects perturbatively. To such an analysis we turn in the 
next section. 



IV. (T-MODEL FOR INTERACTING BOSONS 

To go beyond a non-interacting, semiclassical analysis 
and include the effects of interactions and quantum inter- 

I 



ference, we construct the NLcrM which accounts for fluc- 
tuations of Q. We shift Q by ie„-|-%/2r e"'*^ X e+'^^ + 
^/2Tc e^'-i^ ^{^c + Xj) e+*'i^ to remove these terms 
from the trln[-]. Then, we expand the trln[-] about the 
saddle point and integrate out X, Xc- We obtain an 
effective action which is essentially the same as Finkel- 
stein's action for the fermionic problem (see Appendices 
and I): 



ScffiQ] = J d'^x |L>tr(VQ)^ -4iZtr(eQ) 

ni , . . .,n4 



(12) 



Jaa' Jbb' 



^AA' ^BB' 



P+i-na (^nam 

^aa,BB' ^ 



g i -4 



.BB 



, e '-1' 



where Jab = -^(^ab - a-\g) and S^^ = a\g ± ia\g 
express the particle-hole matrix structure for the density- 
density and Cooper channels, respectively. The parame- 
ter Z is 1 in the bare action above; however, this quantity 
is renormalized, so we have introduced it explicitly here. 
We have absorbed the density-of-states into the diffu- 
sion constant D (and also the coefficients of the other 
terms); the resulting quantity is just the bare conduc- 
tivity and is given by _D = l/27r'' in the above model. 
However, as we noted earlier by considering a model with 
anisotropic masses, and a sufficiently large or small ratio 
mi/ 1712, the bare conductivity will be large. The resis- 
tivity g = l/(27rD) is the expansion parameter used in 
our RG equations, so this observation gives us a limit in 
which they can be applied without apology. 

It may strike the reader as strange that we are using 
a NLcrM to describe a critical point; usually NLcrMs are 
used to describe stable phases because they are so highly 
constrained by symmetry. However, the NLcrM of eq. 
is not, in fact, so rigidly constrained at all. The inter- 
action terms and the tr (eQ) term explicitly breaks the 
0{{k + l)N,kN) 'symmetry' of the model. The latter 
breaks it in such a way as to push the theory into a dif- 
fusive metallic state. However, this symmetry-breaking 
'field' is small in the low-energy limit, so other symmetry- 
breaking fields (or anisotropics) can intervene instead. 
When $ orders in eq. H, Q is forced away from the dif- 
fusive 'direction' in its saddle point manifold, and into 
the superfluid 'plane', where Q has no n- vanishing com- 
ponents which are off-diagonal in particle-hole indices. 
Thus, we can understand the perturbations which lower 
the symmetry of the saddle-point manifold as perturba- 
tions which drive the system away from criticality. There 



are a variety of ways in which one can imagine driving 
the system into an insulating phase. In the absence of 
a better understanding of the Bose glass phase, we con- 
sider the simplest which is just a 'mass' term of the form 
tr {M Q), with M a constant matrix say in replica space, 
which breaks the symmetry of the saddle-point manifold 
and leads to an insulating state. Such a perturbation 
differs only in index strucure with the one imposed by a 
finite Such a term is also generated by shifting /i out 
of the trln[.] term when considering replica symmetry 
broken saddles. Note that none of these possibilities can 
occur in the non-interacting problem, where the symme- 
try of the saddle-point manifold is a genuine symmetry. 

We parametrize Q about the non-interacting saddle 
point as 

Q-'^i'^^^'''^^'^ %,vO (13) 



where the block structure is in frequency space, i.e., the 
matrix is such that n > and m < 0. 

The resulting action is very similar to the 0{N) sigma 
model which is appropriate for a system of fermions with 
spin-orbit scattering. Indeed, one can be transformed 
into the other by redefining q q, q'^ ^q'^ , and 
D —D. The interaction terms look somewhat strange 
at first glance, but the extra i's in (03) are precisely com- 
pensated by the explicit factors of e^^i"* in Eq. (|l2|) (see 
appendix |^) . 
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V. RG EQUATIONS 

Taking advantage of the observation at the end of the 
previous section, we can obtain the RG equations for our 
cr-model by flipping g — > —g in the equations for the 
corresponding fermionic model. Some factors of 2 will be 
different because our bosons are spinless. More details 
may be found in appendix 

The RG equation for is: 



(14) 



Observe that Vc flows to zero, even if 5 = 0. Hence, we 
set Fc to its fixed point value of zero and consider the RG 
equations for g, F, and Z in its absence. To order g^ and 
all orders in F (although, of course, we cannot access non- 
perturbative effects associated with saddle-points which 
are far from the non- interacting diffusive one), the RG 
equations are: 



dg 
M 
dZ 

U 
dr 



^2 2 

= 25 -9 



2 + 2(f-l) ln(l-| 



(15) 

(16) 
(17) 



The physics of these equations is clear from the discus- 
sion in the introduction. Interactions always enhance the 
conductivity to order g^ because the exchange term has 
the same sign as the direct term (they are folded into a 
single F in the bosonic NLaM (^). The gist of the ef- 
fect can be seen from the Hartree and Fock diagrams for 
the boson self-energy displayed in flg. |^. In the Hartree 
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FIG. 1: The Hartree and Fock diagrams for the boson self- 
energy 

diagram, the boson line is repelled by the boson bubble 
which is a measure of the ground state density. (In a 
pure system, this is uniform and cancelled by the neu- 
tralizing background.) In a fermionic system, the Fock 
diagram comes with the opposite sign, so it is an effective 
attraction. In a bosonic system, however, both diagrams 
come with the same sign and lead to a repulsion of par- 
ticles from regions of high density - which, of course, are 
precisely the regions where there are deep wells in the 
random potential. 



The interaction strength, F, grows in importance at 
low energies because it plays a role somewhat analogous 
to the Pauli exclusion principle: in its absence, all of the 
bosons would sit in the lowest minimum of the random 
potential. Z must follow F in order to maintain a finite 
compressibility. 

Notice from Eqs. ( p^ , p^ that Z — T remains invariant 
under the RG flow, as a result of Ward identities that 
originate from charge conservation. It is very useful to 
introduce the coupling constant 7 = F/Z, which allows 
us to rewrite the RG equations in a simpler way: 



dg 

d£ 
dj 

H 



^2 2 

oj -3 



2 + 2- — ^ ln(l-7) 



57 (1 - 7) ■ 



(18) 
(19) 



For (7 > 0, it follows from Eq. (|l8|) that there are two 
flxed-point values 7* = 0, 1 (a closer analysis rules out 
the possibility of another value of 7* with g — 0), as 
shown in Fig. The 7* = flxcd point is unstable, 
while the 7* = 1 one is stable. Consider the RG equa- 
tion for g. The first term on the right-hand-side is the 
weak-localization correction, while the second term is the 
interaction correction. The value 7 = 0.42316 . . . sepa- 
rates the regime where the weak-localization correction 
dominates over the interaction contribution {dg/di < 
for 7 < 0.42316 ... and dg/di > for 7 > 0.42316 . . . ). 
Although the entire surface g = with arbitrary 7 is 
left invariant under the RG flow, any system with bare 
(7, 7 7^ will necessarily flow into the 5 = 0, 7 = 1 fixed 
point. This is the case for short-range interactions, where 
the flow starts with a value 7 < 1. Note that if the bare 
interaction is weak, 7^1, then the resistivity will ini- 
tially increase before eventually decreasing to zero. 

Now, consider the case of dynamically-screened 
Coulomb interactions. As in the fermionic case, the Ward 
identity for charge conservation requires the density- 
density correlation function to vanish at q = 0. This, 
in turn, requires the q-dependent interaction T{q), which 
generalizes F tO[^f^p,case of Coulomb interactions, to sat- 
isfy the identity 



dn 



(20) 



Taking the (7 — > limit o f (pO| ), we obtain Z — T. Substi- 
tuting this identity into (|15|), we see that the second term 
inside the square bracket in ( p^ vanishes. Thus, the RG 
equation for g is dg/di — —3(7^/2, and the resistivity 
flows logarithmically to zero. The system is controlled 
by the same inflnite-conductivity fixed point as in the 
short-ranged case. 

Before concluding this section, let us write down the 
asymptotic behavior near the fixed point g* = 0,7* = 1, 
which we will need later to obtain the critical exponents: 
g - 2/3^ and 1 - 7 ~ exp(- J di g) ^ £"2/3. 
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FIG. 2: RG flow for the resistivity g and interaction parame- 
ter 7 = V/Z 

VI. CRITICAL BEHAVIOR 

The most striking conclusion about the critical behav- 
ior of this system is that the critical resistivity is zero! 
In other words, the 2D superconductor- insulator tran- 
sition is broadly similar to the ZD one. This is some- 
what unexpected. In models such as the Bose-Hubbard 
model, which describes a superfluid insulator transition 
in a clean system, or the 2 -I- 1-dimensional XY model, 
one finds a* — ce^/h, with c a finite universal num- 
ber. At our fixed point, c — oo. Another odd feature is 
the logarithmic approach to the critical resistivity which 
we find; this logarithm is rather different from the type 
which are encountered in the lower critical dimension of 
a phase transition (which happens to be c? = 1 for the 
superfluid-insulator transition). Since a logarithmic flow 
is rather slow, it may not be possible to observe c = go. 
Instead, the critical conductivity at a given temperature 
may actually appear to be a non-universal number which 
depends on the bare conductivity. 

Let us also consider the single-boson density of states, 
N(uj). This may be studied by introducing a source term 
for tr(A(5) into the effective action and computing its 
renormalization. In a system with short-ranged interac- 
tions, we find: 



InN — ~ g ■ "f ■ ln(l — 7) 



(21) 



Substituting the asymptotic forms of g and 7, we find 
that the single-particle density of states diverges weakly, 
A^(<j;) ~ eal'"'"^^/""^-" . Since the boson creation operator 
is the order parameter for the superfluid phase and 



the scaling relation for N{llj) implies that the critical ex- 
ponent 77 = with logarithmic corrections. However, in 
the presence of dynamically-screened Coulomb interac- 
tions, there is a more severe divergence, and we find: 



— In TV = g-e 
d£ ^ 



(23) 



N{uj) = Im (■0t(x, Lv) tp{y., -w)) 



(22) 



Consequently, the single-boson density of states diverges 
at the transition with the power-law N{u!) ~ w"^/'^. This 
implies that the critical exponents 77 and z satisfy r\l z — 
—2/3. Note that we have calculated the density-of-states 
at a metallic critical point. Thus, we should not expect 
Coulomb gap physics to suppress it and give 77 > 0. In the 
fermionic case, the suppression of the density-of-states is 
due to the dominance of the exchange interaction. 

Our NLcrM does not explicitly include single-boson op- 
erators. We assume that their properties can be deduced 
from the the density-of-states. It is certainly possible for 
single-particle operators to be critical even in a theory in 
which only collective modes are retained; this is the idea 
behind bosonization. It is conceivable, however, that our 
NLcrM is incomplete, as regards single-boson properties. 
This could occur if the critical exponent controlling the 
correlation function (7/^^(0;, 0) ?/'(x, 0)) were unrelated to 
that controlling (7/)t(0, r) V'(0, 0)). 

Since Z diverges only logarithmically, the dynamical 
exponent, z = 2, as in a non-interacting system. How- 
ever, in the case of dynamically-screened Coulomb inter- 
actions, there are actually two different diverging time 
scales. One, with exponent z, is the scale associated 
with Z\ it controls the scaling of the specific heat and 
energy diffusion. There is a second exponent, Zc, asso- 
ciated with Z — r, which controls charge diffusion. By 
the same argument as in a fermionic systemtl, eq. |20| 
implies at small q that Z — T ^ q, from which wc con- 
clude that Z — T ^ i.e. Zc = 1. This result was 
obtained for the superconductor-insulator transition by 
a closely-related argument in ref. ^. Combining this with 
our density-of-states calculation, we have rj = —2/3 for 
Coulomb interactions. Notice that ry = —2/3 < sat- 
isfies the lower bound < 2 — d of Ref. O for d = 2. 
The density-of states and the dynamical exponent, Zc, 
are the only quantities which distinguish short-ranged 
and dynamically-screened Coulomb interactions in the 
infrared limit. 

As we discussed in section IV, the leading perturba- 
tion of our a model is a tr (M Q) term, where M is a 
constant matrix say in replica space, which breaks the 
replica symmetry of the diffusive saddle-point manifold 
possibly in the direction of the Bose-glass phase. This is 
a dimension 2 operator at tree level. (If the matrix M is 
proportional to the identity in replica space, this operator 
is instead just a constant at the diffusive saddle-point.) 
Since the coupling constant g fiows to zero, we expect 
a critical exponent = 1/2, up to logarithmic correc- 
tions. This value of v - the mean-field value - violates 
the bound i' > 2/d of ref. 17, However, such violation 



has been seen in other systems as well, and it has been 



8 



are; uedi that the exponent bounded by the theorem of 
ref. |l^ is, in fact, a finite-size scahng exponent which can 
be different from v. 



VII. DISCUSSION 

Diffusion in two dimensions is marginal, and small cor- 
rections (in the limit of large conductivity) such as that 
due to quantum interference or interactions can tip the 
balance one way or the other. Contrary to conventional 
wisdom, it is hardly a foregone conclusion which effect 
will win. After all, weak localization is weak. Interactions 
can easily overpower it, leading to metallic behavior. Ac- 
cording to our analysis, this is precisely what occurs at 
the superconductor-insulator transition. The effect of in- 
teractions is so dominant that the universal value of the 
conductivity at the transition is infinity. Such a diverging 
conductivity has been found in models, with interaction 
and dissipation, but without disorderO 

The possibility of a metallic phas|e-|j™thin the Bose 
glass phase has been studied recentlyc3l2il. We focus on 
the diffusive properties at the critical point, and do not 
investigate whether saddle-point solutions whithin the 
Bose glass could lead to non-zero conductivities. How- 
ever, it is noteworthy that an infinite critical conductivity 
is consistent with a Bose metal with a diverging conduc- 
tivity at the transitioncil. 

We derive these results in a NLcrM approach, in which 
we discard those critical modes of the clean system which 
are extraneous and retain only the particle diffusion 
modes of the disordered system. The resulting NLctM 
leads to a number of non-trivial predictions: (1) the crit- 
ical conductivity is infinite; (2) there are two diverging 
times scales if the interaction is Coulombic, one associ- 
ated with charge diffusion, which has exponent 2 = 1, the 
other associated with energy diffusion, which has expo- 
nent z = 2; (3) the single-boson density of states diverges 
as uj^"^/^ , which implies a critical exponent rj — —2/3 in 
the case of Coulomb interactions; for short-range inter- 
actions, it diverges logarthmically; (4) the correlation- 
length exponent takes the mean- field value v — 1/2. 

If boson-vortex duality were to hold exactly, then one 
would expect <?* = 1 (in units of (2e)^//i). Our result ap- 
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pears to imply that duality is violated logarithmically: 
bosons are more mobile than vortices in the infrared 
limit. However, it is hard to see how the physics of vor- 
tices enters at all into our calculation, so it is possible that 
we have missed important non-perturbative effects. Our 
results, do not agree with the numerical study of Wallin, 
et aZ.E3. However, the flow to our fixed point is logarith- 
mic, and this may be too slow for a numerical study on 
a finite-sized system. Alternatively, they may simply be 
accessing a different fixed point which attracts systems 
with small bare conductivities. And finally, since their 
starting point studies phase but no amplitude fluctua- 
tions, the two models may simply be in different univer- 
sality classes. Our results also differ quantitatively from 
those Herbut, which are based on an expansion about 

The measured critical exponents for the zero-field 
superconductor-insulator transitio n , .. wh ich is accessed by 
varying the thickness of a thin filmc3l2a, are those of clas- 
sical percolation. This does not agree with our theory, 
but it also suggests that the experiments are not quite 
in the asymptotic quantum critical regime, but rather in 
some higher-temperature classical regime. There is dis- 
agreement about the values of the critical exponents at 
the magnetic-field-tuned rSuperconductor- insulator tran- 
sition. One experiaaentEa finds percolation-like expo- 
nents, while anotherCJ finds v — 0.7 ±0.2, which includes 
our theoretical prediction at the edge of its error bar. 
(All of thesCptprperimcnts find z 1, as expected on gen- 
eral groundffla, and in our theory.) The applicability of 
our strategy to a magnetic-field-tuned superconductor- 
insulator transition is a question for future study. 



Acknowledgments 

We would like to thank M.P.A. Fisher, V. Gurarie, 
A.W.W. Ludwig, and P. Phillips for discussions. C.N. 
was supported by the National Science Foundation under 
Grant No. DMR-9983544 and by the A.P. Sloan Founda- 
tion. C.C. was supported by the National Science Foun- 
dation under Grant No. DMR-9876208 and by the A.P. 
Sloan Foundation. 



APPENDIX A: DERIVATION OF THE cr-MODEL 

Here we derive the effective action for the interacting disordered bosons in terms of two fields $ and Q. We work, 
in sequence, on the free part, the disorder part, and finally the interaction part of Eq. (^. 



The free action 



We start by introducing a bosonic amplitude $ to decouple the chemical potential (/^) term. $ acquires a finite 
expectation value when bosons condense. 
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The free part of the action 



S'frcci^] = ^ J (fxi (j)naAix) Ne„ + + ^ j A„™ 4>m.aA{x) (Al) 

can is generated upon integration of a decouphng field $„(i = $fiai + 'i-^na2 in 



Sfrcoi^, <&] f (fx i (jjnaAix) Ne,i + ) 



0mayl(a;) 



2. Disorder term 

Let us next decouple the four bosons in the disorder term in Eq. 

•Srand = ^ J (f X ^ 'PnaA{x) Ann' (pn' aA{x) (j),nbB{x)A,nm' 4>m' bB{x) (A3) 
n,n' ,rn^m' 

where Amm' = sgn(m) (5„im' ■ The same disorder term is generated upon integration of the Hubbard-Stratonovich 
matrix field Q™^^ 

= JdQ e-^ f ""^ e-^«^[«'^] (A4) 

S'HsiQ, </>] = « X! J d^X (j)naA{x) QaJl^ABix) A„i„i>(j)m>bB{x) (A5) 



where 



The matrix Q has indices in three separate spaces, i.e., it is assembled as a direct product in energy n, m, replica a, 6, 
and real-imaginary A, B spaces. The trace of corresponds to 

= QTab QZ-ma , (A6) 

where repeated index summation is carried out in all three spaces. When we write for short Qnn' we mean a matrix 
whose elements are matrices in replica, and real-imaginary spaces. 

3. Interaction term 

Let us consider the case of short range interactions, in the density-density (s) and pairing (c) channel. Once again, 
we will omit sums over indices for replica and real-imaginary parts, and write explicitly the Matsubara sums. 

5i„t = Ss + S, (A7) 

where 



Ss = Ts ^ d'^xe [4>niaA{x) JaA' 4>n2aA{x)][<j)n3aB{x) JbB' (pmaBix)] Sni-n2+n3-n4 (A8) 

5'c = Lc J d'^Xe^'T^^''^'^^"^^ [(j)n,aA{x) S^^, (t)n2aA'ix)][(j)n3aB{x) Sgg, (f>n4aB'{x)] <5„i+„,_„3_„, (A9) 

ni ,. . . ,714 

with the matrices Jab = -^{^ab — <^ab) ^^'^ ^ab — '^ab ^ ^'^\b i^^^ being Paufi matrices). Notice that the 
different terms within square brackets above correspond, in terms of the original bosons t/'j to E^iid V'V'- 
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We now introduce two Hubbard-Stratonovich fields X and Xc,Xc* to decouple the four <j) interactions: 



(AlO) 



^HS/xW = J2 J ct>naA{x) e-^f^8"(") XXab(^) e"^* ^^'^("VmbB (x) (All) 



and 



DXr e"-^^!^"! e""^"s/Xe[^o 



e 4 



(A12) 

^*^"('"Vm6B(2;) (A13) 



where 



^ab,AB — -^a "ab -JAB ^Cab,AB — "ab '^AB 

Notice that the matrices X"^^^ and Xc^b^y^g depend, respectively, only on the energy difference n — m and sum 
n + m. The action for the matrices X and Xc is 



SAX] = l^fd^xx: 



(A14) 
(A15) 



4. Integrating out the <p fields 



We can summarize all terms above: 



X,= J Dcj) D^DQDXDX,*DX, e"^ ^ '^'^ e"^-!^! e-^-o[Xo] 

^g-5free[0,*] g-SHs[0,Q] g-SHS/x[0,X] g - SrS/X 

where we can express the *S'free5 "Shs, S'hs/x, and S'hs/x„ in a more concise (matrix) notation as follows: 
&s[0,(3] = j d^x f [x) iQ{x) K(t>{x) 



'HS/X 



/2I\ j d^x (jF{x) e-'^i^ X{x) 



-if A 



<i>{x) 



HS/Xo 

where the matrix = e„ (5nm- 

Integrating out the boson fields xj), we obtain 



[X,] = iV2rr d'^x(t)^{x)e-''i^ - X^{x)+xXx) e-^^A^^^.) 



with 



X, = Jd^P D^DQDXDX,*DX, 1 '^'^ e"^--!^! e-^-[^=l e-^°[«'*'^'^= 

S'o[Q, X, = II J d^x ^'^ G ^ + J d^x tr log i (^ifi + ^V^^ A + iQ{x) A 

+i^/2Fse~'^^ X{x) e-'^^ + i^/2F^e-'^^ ^(^Xa{x)+Xa\x)^ e'^^' 



(A16) 

(A17) 
(A18) 
(A19) 

(A20) 

(A21) 



(A22) 



iJd'^x^^G^ + J d'^x trlog (^m + ^V^^ + Q{x) 
VstT e-'^^ X{x) e+'^^ - e''^^ i (^X^ix) + Xj{x)^ e+'^^j + const. (A23) 
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The propagator G depends on Q, X, and Xc, X*. 



5. Shifting Q 

Let us now shift Q: 



i^/2Tr, e-*TA z.[x^ + X^ 



and 



trg2 = tr Q2 



+i2j2Ts tr 



+ i2V2rc tr 



(A24) 
(A25) 



The Matsubara cut-ofF k comes from the extra frequency sum in the trace, and there are factors of 2 from the traces 
over the real-imaginary components, tr J = 2 and tr S~^S~ = 4. 

The next step is to integrate out the X, Xc fields. This generates quadratic 

in Q terms (we now drop the tildes for notational simplicity). It is useful to define 



" ~ ,,2 



Vo 1-8 



Iaa'bb' 



„+«f sgn ria /-)n3«4 -if sgn n4 
C '^aa.BB' ^ 



where the tensors 7aa',sb' depend on the channel: 

7aa',bs' = Jaa' Jbb' 

i'aa',bb' = Sgg, 

Summarizing it all, we have an effective action 



S'e//[Q,$] = ^ j d'^xtiQ^ + j d'^x trlog (^in + ^V^^ + Q{x) + ^ J d'^x <^'^ G ^ 



^aa,AA' ^ 
g+if sgn m r)"l"2 g-if sgn 712 



^+»f sgn TI3 --if sgn ^4 

^ ^aa,BB' ^ 

^+»f sgn ns /^n3"4 -if sgn 714 
6 ^aa,BB' ^ 



(A26) 

(A27) 
(A28) 

(A29) 

^ni— 712+^3— 714 



^711 



+712—713—714 



r 



APPENDIX B: SEPARATION OF 
LONGITUDINAL AND TRANSVERSE 
(DIFFUSIVE) MODES 



In this appendix wc show that the transverse fluctu- 
ations of the Q field correspond to boson diffusion for 
H = 0, similarly to the fermionic case. We show that 
the diffusion term arises even in the absence of a small 
parameter 1/EpT. 



Expansion to quadratic order in 5Q leads to a term 
\ M-^-^iq) SQZ^^A-Q) ^Q7::a'a{q) (B1) 



where 



_ 1 _ f d'^p 
~Vq~ J {2nY 



G{p,ni)G{p + q,n2) . (B2) 



The first term on the right-hand side comes from the 
trQ^ in the action, and the second term has its origin 
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in the trlog(-). M"^'^'^{q) is selected to be diagonal in 
and independent of replica a, a' and real-imaginary A, A' 
indices. The Green's function 



G{p,ni) 



1 



ie„ - E{p) + 27 sgn(e„) 



(B3) 



For sgnm sgnn2 > 0, the real part of the integral vanishes 
for (? ^ 0, so that Re[ M"i"2(g)] = ^, and we are left 
with a massive longitudinal mode. 

Let us turn to the interesting case sgnni sgnn2 < 0. 
For simplicity, we neglect the itni 2 terms in the denomi- 
nator (these terms can be handled alternatively by shift- 
ing the Q field). Expanding the integral in Eq. (B2) in 
powers of q: 



dfip 



df^p 



(27r) 



G+{p) G-ip) 



1«: 

2 m 



d'^p 



Glip) G-(p) 



4Eip) 



GUp) G_(p) 



(B4) 



where G±{p) = [—E{p) ± ^]^^ ■ The integrals over mo- 
menta can be transformed into integrals over energy e 
using the density of states v{() — (2m)''/^e''/^~^ 
Define 



Ia,b.c= r A 1/(6)6'^ [G+{et [G^{e)r 
Jo 



(B5) 



so 



(2^) 



-^G'>ip,m)G'ip + q,n2) 



I +^'^' 



0,2,1 



4. 

3^1,3,1 

a 



(B6) 



(a finite upper frequency cut-off is needed depending 
on d, a, b, c). It is also convenient to rescale the energies, 
defining y = 2Te, so we can write 



2 (2n)° 



X 



2rf2 



dy 



d/2-l+a 



Ad (2t) 
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-a-\-b-\-c 



i-y + i)'' i-y-iY 



2rO 



dy 



y 



d/2-l+a 



(B7) 



with Ad 



1 



2 (27r) 



- (2to)'^/2 (2r)-'*/2 



One can check that once r is fixed by the saddle point 
Eq.(no[), which can be cast as 



Im /( 



0,1,0 



-Im /( 



1 



0,0,1 



2tvo 



(B8) 



then it follows trivially that 

-^0,1,1 = 



1 



(B9) 



so that the leading order term in AI±{q) is of order q^, 
which allows us to define the diffusion constant 



D 



1 

Am 



Io,2,i + -ih.s,! 
a 



(BIO) 



The last step remaining is to show that D is purely real. 
After simple manipulations, one can show that 



ImL> = 



J Sd_ 

4m (2n)<^ 

X / dy 





{2mf^^ (2t)- 



-d/2+3 



1 2/ - 1 



y 



d/2-l 



(y2 + 1)3 



(Bll) 



It is trivial to show by integration by parts (splitting 
the integrands into f{y) = y/{y^ + 1)^ and g{y) = y") 



that the integral in Eq. ( pill ) scales as (t17)'*/2"^. Thus 
the cut-off can be safely taken to infinite for d < 8, and 
Iml? = 0. 

Notice the difference between the fermionic and 
bosonic cases. In the fermionic case one can also in- 
terchange momentum p integrals for energy e integrals, 
using the density of states at the Fermi level Ep- The in- 
tegrals are cut-off by the bottom of the band, —Ep away 
from the zero energy states. In the bosonic case, one 
starts from the bottom of the band, and needs to include 
an energy dependent density of states i^(e); the cut-off H. 
is introduced only for convergence, and 51 — > cx3 is possi- 
ble for d < 8. In contrast to the fermionic case, where Ep 
is finite, in the bosonic case for a perfect parabolic spec- 
trum Q oo. The small parameter for the Fermi case is 
{Eft)"^, whereas for the Bose case it is {ftT)~^ — > 0. 



APPENDIX C: PARAMETERIZATION OF THE 
SADDLE AND RELATION TO THE FERMIONIC 
cr-MODEL 

As we previously mentioned, we can easily obtain 
the RG equations for the conductance and interaction 
couplings by determining a correspondence with the 
fermionic model. Here we show how this is achieved. 

Let us first look at the Finkclstein type terms in the ef- 
fective action for the Q fields. The Q matrices are param- 
eterized as in in Eq.(|l3|), repeated here for convenience. 



Q 



mvo 




qq^Y 



-i{\ + q^q) 



(CI) 



The quantities that appear in the Finkelstein type terms 
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for the bosonic problem are 

mvo I i (l + qq^^ ^ iq 

-i{q^) -I (1 + q^q) 

{\-q{-qT)Y q 




1 




Direct comparison with the fermionic saddle point 



[l-qq^Y q 



q 



(C4) 



shows that the terms in the Finkelstein type action for 
bosons, upon parameterization in terms of q, q'^ , are the 
same as the ones for fermions upon the identification q 
q and q'^ — > —< 



The extra factor of 



03 



squared (because the Finkelstein terms are quadratic in 



Q) , makes the sign of the interaction term for bosons and 
fermions the same. 

For discussing the diffusive term D j d'^xtclVQY , no- 
tice that by rewriting 



2 \^{-qr) 



. (C5) 



we again identify it with the fermionic saddle point Qf^ 
but now the off-diagonal elements have extra factors i, —i. 
These factors will cancel each other in the expansion 
of the quadratic in Q diffusive term, and hence can be 
dropped, and once again the fermionic saddle expansion 
can be used. The overall factor of i has the effect of 
changing D ^ —D. 

In summary, all the RG equations for the dirty inter- 
acting boson problem can obtained from those of the (in- 
teracting) fermionic orthogonal ensemble upon replacing 
g —g (or D —D). 
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